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Abstract 


For  conventional  imaging  systems,  geosynchronous  earth  orbit  (GEO)  space 
objects  cannot  be  resolved  due  to  their  40  Mm  distance.  There  exists  a  strong  need  to 
obtain  high  resolution  images  of  GEO  objects  and  to  accomplish  this  task,  investigation 
into  the  suitability  of  inverse  synthetic  aperture  laser  radar  (ISAL)  is  currently  underway. 
A  critical  component  in  detennining  this  suitability  is  to  accurately  model  the 
atmospheric  impacts  on  laser  detection  and  ranging  (LADAR)  pulses.  The  expected 
result  is  that  atmospheric  impacts  are  a  strong  function  of  the  illumination  beam  length. 

Conventional  knowledge  says  that  while  the  atmosphere  chums,  wind  is  the 
predominant  cause  of  temporal  evolution  which  simplifies  all  modeling  and  simulation 
into  the  frozen  flow  hypothesis.  The  concern  is  that  the  frozen  flow  hypothesis  based 
phase  screen  generation  techniques  fail  to  accurately  predict  the  temporal  development  of 
optical  phase.  Additionally,  wavefront  phase  aberrations  introduced  by  the  atmosphere 
are  dominated  by  tilt,  which  is  often  underrepresented  in  the  frozen  flow  approximation 
and  causes  the  most  degradation  of  direct-detection  LADAR  system  performance. 

This  thesis  proposes  a  new  approach  and  provides  a  detailed  derivation  of  a  new 
temporally  evolving  Zernike  polynomial  based  atmospheric  phase  screen  generation 
model.  This  new  model  is  verified  through  laboratory  experiment,  and  then  utilized  to 
analyze  atmospheric  effects  on  long  pulse  LADAR  beams  and  their  impact  on  mixing 
efficiencies.  It  is  shown  that  this  new  turbulent  flow  model  more  accurately  predicts 
mixing  efficiency  than  that  of  the  basic  frozen  flow  approximation. 
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Atmospheric  Impact  on  Long  Pulse 
Laser  Detection  and  Ranging 
(LADAR)  Systems 

I.  Introduction 

Imaging  otherwise  non-resolvable  Geosynchronous  Earth  Orbiting  (GEO)  objects 
with  Inverse  Synthetic  Aperture  Laser  detection  and  ranging  (LADAR),  or  ISAL  is 
currently  the  focus  of  a  recent  Air  Force  Research  Laboratory  (AFRL)  effort.  ISAL 
acquires  images  of  a  coherently  illuminated  target  as  it  traverses  the  sky.  This 
synthesizes  a  much  larger  aperture  than  that  of  the  observatory  and  through  post 
processing,  permits  the  construction  of  high  resolution  images.  Considering  the  3.67  m 
primary  mirror  of  the  Advanced  Electro  Optical  System  (AEOS),  and  a  typical  operating 
wavelength  k  of  0.5  jum,  the  minimal  resolvable  angle  Ae  is  equal  to  Ae  =  1.22A/3.67, 
or  roughly  0. 1667  //rad.  Transposing  that  to  the  smallest  resolvable  object  at  a  typical 
GEO  distance  of  40  Mm  equates  to  approximately  6.5  meters.  The  largest  of  GEO 
satellites  therefore  occupy  no  more  than  a  few  pixels  of  an  image.  ISAL  is  intended  to 
overcome  the  physical  limitations  of  meter  class  telescopes  such  as  the  AEOS. 

1.1  Problem  Statement  and  Hypothesis 

This  thesis  evaluates  the  atmospheric  impacts  on  long  pulse  LADAR  systems,  a 
manageable  portion  of  the  much  larger  ISAL  problem.  More  specifically,  this  thesis 
analyzes  mixing  efficiencies  as  a  function  of  the  pulse  length.  The  analysis  provided 
herein  aids  in  the  investigation  of  ISAL. 
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Conventional  knowledge  implies  atmospheric  impacts  are  a  strong  function  of  the 
illumination  wavelength  and  beam  duration.  The  often  ignored  additional  dimension  of 
temporal  evolution,  the  boiling  and  churning  of  turbulence,  is  the  primary  focus  of  this 
thesis.  Even  during  the  relatively  long  beam  pulses  used  in  ISAL,  the  atmosphere  simply 
does  not  have  time  to  “shift”  according  to  the  frozen  flow  hypothesis  approximation. 
Note  that  ISAL  utilizes  a  chirped  waveform,  which  is  paraphrased  as  a  pulse  within  this 
thesis.  Consider  for  a  moment  a  100  //sec  long  pulse  where  the  combined  and  prevailing 
wind  is  10  m/sec.  If  prognostic  analysis  is  limited  to  frozen  flow,  then  the  atmosphere 
will  shift  a  single  millimeter  during  that  pulse.  For  shorter  pulse  lengths,  lower  wind 
speeds,  or  a  combination  of  the  two,  the  superficial  shift  is  even  less  significant.  It  is 
expected  that  the  affects  of  the  atmosphere  upon  the  pulse  are  directly  related  to  the 
illumination  wavelength  and  the  duration  of  the  pulse. 

1.2  Research  Focus 

It  is  hypothesized  that  since  the  atmosphere  cannot  flow  by  any  considerable 
amount  on  these  time  scales,  perhaps  it  can  boil,  or  chum.  Therefore,  the  critical 
component  of  this  thesis  is  the  development  of  a  new,  accurate  atmospheric  phase  screen 
generation  technique  that  not  only  shifts  as  the  frozen  flow  model  does,  but  also  evolves. 
This  new  model  permits  realistic  atmospheric  evolution  during  the  relatively  short  pulses 
of  spotlight-mode  ISAL  illumination,  and  ultimately  the  detection  and  imaging  of  GEO 
space  objects. 

The  predominate  focus  of  this  research  is  the  development  and  validation  of  the 
new  temporal  phase  screen  generation  model.  As  previously  mentioned,  it  is  expected 
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that  atmospheric  impacts  are  a  function  of  the  ISAL  system’s  pulse  lengths.  This  new, 
unique  and  useful  tool  should  not  only  verify  such  conventional  understanding,  but  can 
be  used  in  numerous  applications,  particularly  within  LADAR  based  applications.  This 
model  is  validated  through  laboratory  experiment  in  which  recorded  data  will  be 
collected,  then  compared  to  theoretical  as  well  as  the  new  and  commonly  used  models. 
The  modeling  and  simulation,  recorded  lab  data  analysis,  and  subsequent  validation  are 
perfonned  within  the  Mat  lab3'  environment. 

1.3  Thesis  Overview 

Following  this  brief  introductory  chapter,  Ch.  II  provides  an  introduction  to  ISAL, 
simulating  atmospheric  turbulence  with  phase  screens,  and  mixing  efficiency.  Two 
principal  existing  methods  of  phase  screen  simulation  are  discussed:  Fourier-based  and 
Zernike-based  approaches.  The  development  of  temporal  correlation  of  the  new  Zernike- 
based  turbulent-flow  phase  screen  generation  is  the  topic  of  Ch.  III.  The  laboratory 
configuration  which  provides  measured  data  for  the  new  modeFs  comparison  is  discussed 
in  Ch.  IV  as  well  as  the  results  of  this  with  comparative  analysis.  The  methodologies  of 
simulation  used  to  achieve  the  objectives  of  Sec.  1.2  are  discussed  in  Ch.  V,  as  well  as  a 
presentation  and  discussion  of  the  findings.  Finally,  Ch.  VI  summarizes  the  research 
efforts  of  this  thesis,  discusses  the  challenges,  critical  results,  and  provides  guidance  for 
future  related  research  efforts. 
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II.  Background  Theory 


This  chapter  introduces  the  basic  concepts  and  components  of  the  ISAL  system 
that  are  modeled  into  simulation  for  subsequent  analysis.  Basic  phase  screen  generation 
techniques  are  discussed  in  terms  of  simulating  atmospheric  turbulence,  with  particular 
emphasis  on  the  Zemike  polynomial-based  approach.  This  is  because  Ch.  Ill  develops 
the  new  model  based  on  the  basic  Zemike  theory  provided  here.  Additionally,  an 
overview  of  the  system  is  presented,  highlighting  the  relevant  concepts  of  interferometry 
and  its  application  to  mixing  efficiency. 

2.1  Inverse  Synthetic  Aperture  LADAR  (ISAL) 

Conventional  imaging  systems  require  a  10m  class  or  larger  aperture  to  provide 
useful  infonnation  of  distant  GEO  objects  [12].  Synthetic  aperture  systems  rely  on 
relative  movement  between  the  target  and  observation  point  to  effectively  synthesize  a 
larger  aperture.  This  permits  cost  savings  by  utilizing  smaller  optics.  Where  synthetic 
aperture  radar  (SAR)  systems  rely  on  a  moving  observer  with  a  static  target,  ISAL 
conversely  relies  on  a  moving  target  with  static  observation. 

ISAL  has  three  key  characteristics  which  differentiate  it  from  these  previously 
proposed  solutions:  tomographic  imaging  techniques,  coherent  illumination  and 
heterodyne  detection,  and  utilizing  a  chirped  ISAL  waveform  [12].  These  three 
characteristics  in  combination  are  expected  to  permit  high  resolution  imaging.  Utilizing 
coherent  illumination  allows  for  tomographic  reconstruction  from  the  Fourier  data,  and  in 
combination  with  heterodyne  detection,  provides  methods  for  boosting  the  signal  to  noise 
ratio  (SNR). 
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2.2  System  Overview 

While  utilizing  ISAL  eliminates  the  difficulty  and  cost  of  obtaining  a  10m  class 
telescope,  it  requires  actively  illuminating  GEO  objects  and  coherently  detecting  the 
return  signal.  The  process  of  active  illumination  incurs  atmospheric  perturbations  of  both 
the  outbound  and  return  waveform.  However  the  atmosphere  effectively  ends  at  600km, 
1.5%  of  the  40  Mm  distance.  As  such,  regardless  of  any  irradiance  fluctuations,  beam 
spreading,  or  loss  of  spatial  coherence  of  the  outbound  wavefront  the  remaining  39.4  Mm 
vacuous  propagation  ensures  that  even  the  largest  GEO  satellites  are  uniformly  and 
coherently  illuminated  [1].  Therefore,  as  the  goal  is  to  determine  the  mixing  efficiencies 
as  the  pulse  length  varies,  the  simulation  may  be  modeled  as  a  one-directional, 
illuminated-GEO-object-to-observation-point  propagation.  Pulse  length  dictates  the 
duration  in  which  light  propagates  through  the  atmosphere,  directly  limiting  the  time  for 
the  atmosphere  to  evolve.  It  is  during  this  period  that  relative  phase  differences  between 
the  transmitted  and  reference  beams  affect  mixing  efficiencies. 

It  is  readily  understood  that  even  the  most  impressive  adaptive  optics  (AO) 
systems,  such  as  those  within  the  AEOS,  cannot  continuously  and  fully  correct  wavefront 
perturbations.  If  the  AO  system  could  fully  compensate  then  there  would  be  no  losses  in 
mixing  efficiency.  This  uncorrectable  time  interval  in  which  modulation  efficiency 
losses  occur  is  what  concerns  AFRL.  Recall  that  if  temporal  evolution  is  limited  to  a 
frozen  flow  approximation,  the  AO  system  will  require  minimal  to  no  compensation  due 
to  nearly  non-existent  shift. 

The  ISAL  system  is  expected  to  be  integrated  into  the  3.67m  class  AEOS  at  the 
Air  Force  Maui  Optical  and  Supercomputing  (AMOS)  observatory,  where  some  of  the 
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best  seeing  in  the  world  is  common.  A  seeing  parameter  (Fried's  constant),  or  ro,  of  14 
cm  is  typical  and  is  the  value  utilized  in  simulating  the  objectives  of  Sec.  1 .2.  The  GEO 
object  will  effectively  be  modeled  as  a  point  source  illuminated  by  varying  pulse  lengths 
with  a  visible  wavelength  X  =  550  mn.  The  GEO  object  illumination  propagates  through 
a  discrete  number  of  independently  evolving  phase  screens  to  the  3.67m  aperture  as 
described  in  Sec  5.1.1.  Then  after  propagating  to  the  detector  plane,  the  measured  test 
beam  is  combined  with  the  control  beam  to  measure  mixing  efficiency. 

2.3  Atmospheric  Model 

Turbulence  is  largely  located  within  relatively  few  narrow  concentrated  layers  at 
good  astronomical  sites  such  as  AMOS  [1].  Consider  for  a  moment  a  critical  component 
of  optical  imaging  systems:  lenses.  While  a  lens  can  converge,  diverge,  focus  and 
disperse  light,  the  most  important  property  of  a  lens  is  the  phase  transfonnation 
introduced  to  the  propagating  light  [6].  A  lens  in  which  light  enters  and  exits  at 
approximately  the  same  transaxial  coordinates  is  referred  to  as  a  thin  lens.  Thin  lenses 
primarily  delay  the  incident  wavefront  proportional  to  the  thickness  of  the  lens  at  each 
coordinate.  In  this  sense,  thin  lens  approximations  of  atmospheric  turbulence  and  phase 
screens  are  synonymous.  This  discrete  modeling  of  the  atmosphere  is  pennitted  by  the 
profile  of  the  structure  constant  of  (h),  the  refraction  index  fluctuations,  with  discrete 
values  of  height  h  [1],  The  layers  of  atmospheric  turbulence  are  simulated  with  random 
phase  screens  whose  power  spectrum  follows  a  Kolmogorov  model. 

Two  common  phase  screen  generation  models  are  comparatively  analyzed  for  the 
objectives  outlined  in  Sec  1.2.:  the  fast-Fourier-transform  (FFT)-based  method  and  the 
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Zernike  polynomial-based  method.  Other  methods  such  as  the  Fourier  series  model  will 
not  be  analyzed  as  they  use  the  same  frozen  flow  temporal  evolution  that  the  FFT  method 
does  [15].  The  FFT-based  method  will  provide  comparison  of  frozen  flow  approximation 
to  the  new  Zemike-based  turbulent  flow  results.  The  basic  Zernike-based  theory 
presented  here  is  expanded  to  provide  temporal  evolution  in  Ch.  III. 

2.3.1  FFT-based  Phase  Screen  Generation 

The  FFT-based  methods  are  quite  common  since  very  large  phase  screens  may  be 
generated  quickly.  FFT-based  methods  are  similar  in  that  they  are  generated  from  their 
respectively  named  power  spectral  densities  (PSDs).  Put  simply,  these  screens  are 
generated  by  taking  the  inverse  FFT  of  the  square  root  of  the  PSD  multiplied  by  a  random 
zero  mean,  unit  variance  matrix.  See  Fig.  1  for  a  basic  visualization  of  this  process. 


y1 


Zero  /i, 
unita, 
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Figure  1:  FFT-based  phase  screen  generation. 


One  well  known  disadvantage  of  this  approach  is  that  lower-order  aberrations  such  as  tilt 
are  often  under-represented  [14].  These  lower-order  aberrations  make  up  a  majority  of 
the  atmospheric  energy  spectrum  and  as  such,  must  be  included  to  produce  realistic 
models.  One  of  the  more  widely  used  FFT-based  methods  is  that  of  the  modified  von- 
Karman  spectrum, 

/  \  2 

=  0.023(r0)-5/3e  ^ Km '  (k2  +  k02)~11//6.  (l) 


The  modified  von-Karman,  allows  an  application  of  the  inner  do)  and  outer  ( Lq )  scales  of 
turbulent  eddies,  leading  to  Km  =  5.92 Ho,  and  kq  =  IkILq,  in  which  k  is  the  spatial 
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frequency  and  is  scaled  by  a  factor  of  Fried’s  constant,  tq.  This  investigation  assumes 
ideal  limits  where  Lo  approaches  infinity  and  lo  approaches  zero,  as  this  best  reflects  the 
inner  and  outer  scales  of  the  subsequently  discussed  Zernike-based  method.  With  these 
inner  and  outer  scales,  Eq.  (1)  reduces  to  Eq.  (2). 

0>(#c)  =  0.023(r0)"5/3(fc2)-11/6  (2) 

A  sample  screen  generated  by  this  method  is  shown  in  Fig.  2.  Discussion  of  generating 
FFT-based  screens  is  limited  to  this  point  as  they  are  well  known  [9],  [14],  [15],  and  only 
utilized  for  comparative  analysis.  Further  discussion  and  methodology  for  generating 
FFT-based  frozen  flow  approximation  phase  screens  is  readily  available. 


Figure  2:  FFT-based  phase  screen. 


2.3.2  Zernike-based  Phase  Screen  Generation 

Drawing  a  stark  contrast  to  Fourier  methods,  the  Zemike  polynomial  approach 
does  not  begin  with  a  random  phase  array.  Instead  coefficients  are  combined  into  two- 
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dimensional  random  functions.  In  his  1976  paper,  Noll  [10]  defined  a  new  modified  set 
of  Zemike  polynomials,  Zj(r,  0).  He  defined  these  polynomials  to  be 

Zevenj(r>  8)  =  Vn  +  1S™V2  cos(m0)) 

Zodd  j  0, 9 )  =  Vn  +  1S™V2  sin(m0)  )  (3) 

Z;  (r,  0)  =  Vn  +  15°  (r),  m  =  0 

where  the  functions  S™(r)  are  referred  to  as  radial  functions  and  are  defined  as 

n-m 

SmM  =  V  _ (-1  H»-s)! _  (4) 

Z_i  s!  [(n  +  m)/2  —  s]!  [(n  —  m)/2  —  s]! 

s=o 

This  function  is  based  on  a  polar  coordinate  system  where  S  is  the  radius,  r  is  the  position 
along  the  radius,  and  6  is  the  angle  with  respect  to  the  x-axis.  Noll’s  numbering  sequence 
of  the  index  j  corresponding  to  each  polynomial  Z;  proceeds  by  row,  then  value  of  radial 
degree  n,  and  increases  with  azimuthal  frequency  m.  Combining  these  equations,  it  is 

then  possible  to  generate  independent,  random  phase  screens  &atm  from  the  relationship 

J 

®atm {Sr,  0)  =  ^  djZj  (r,  0) .  (5) 

;=i 

To  generate  the  phase  screen  with  accurate  correlation  it  is  necessary  to  calculate  the 
coefficient  amplitude  vector  a;  by  multiplying  <f>  with  a  zero-mean,  unit-variance  random 
vector  n.  The  n  vector  has  a  length  equal  to  J.  The  a;-  vector  is  defined  as 

CLj  =  On,  (6) 

where  O  is  the  Cholesky  decomposition  of  the  covariance  matrix, 

0  =  Vc  and  OtO  =  C.  (7) 
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The  covariance  matrix  C  is  generated  from  the  covariance  of  any  two  Zemike  polynomial 


coefficients  a;-  and  a-'  [13]  as  provided  in  Eq.  (8): 


,  r  1  =  _ Kzz-Szr [(n  +  n-  5/3)/2 ](P/r0)5/3 _ 

■j,j'  I  ai’  aj'\  r[(n  _  n'  _  i7/3)/2]f[(n  -  n  -  17/3)/2]f[(n  +  n'~  23/3)/2]' 


(8) 


in  which 


K'  = 


2712 


(n+n  -2m)  _ 

x  (—1)  2  -\/  (n  +  l)(n  +  1) 


(9) 


and 


8Z  —  (m  =  m  )  | parity (J,j ')  V  (m  =  0)].  (10) 

An  example  of  a  phase  screen  generated  by  this  method  using  1 ,024  Zernike  polynomials 
is  provided  in  Fig.  3. 


Figure  3:  Zernike-based  phase  screen. 
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2.4  Mixing  Efficiency 

Mixing  efficiency,  referenced  interchangeably  with  modulation  efficiency  throughout  this 
text,  is  the  percent  efficiency  of  combining  a  test  signal  with  a  control  signal.  It  should 
be  accentuated  that  analysis  of  mixing  efficiency  has  not  previously  been  accomplished 
in  terms  of  IS  AL  application.  The  illumination  signal  is  split  and  then  later  combined  by 
interferometers,  or  beam  splitters,  as  Fig.  4  depicts. 


Figure  4:  Interferometer  functionality. 


The  intended  application  in  regards  to  ISAL  involves  boosting  the  target  return  signal 
SNR  by  adding  it  to  an  undistorted  and  delayed  control  signal.  In  other  words,  mixing 
efficiency  is  the  percent  efficiency  in  which  a  split  beam  propagated  through  different 
paths,  may  be  recombined.  Conventional  knowledge  assumes  that  obtainable  modulation 
efficiencies  are  approximately  10%.  Sec  5.2  provides  the  mathematical  methodology  for 
accomplishing  this  analysis.  For  the  purposes  of  this  introduction,  mixing  efficiency  is 
defined  as  the  fraction  of  test  beam  intensity  minus  the  bias,  all  of  which  is  divided  by  the 
bias;  where  the  bias  is  half  of  the  control  beam  intensity.  The  control  intensity  is 
calculated  by  combining  the  intensity  of  two  vacuous  far-field  propagated  point  sources, 
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subsequently  synonymous  with  control  beams.  The  test  intensity  is  calculated  by 
combining  the  intensity  of  a  vacuous  far- li eld  propagated  point  source  with  the  intensity 
of  point  source  propagated  through  turbulence  [8],  The  intensity  ID  incident  on  the 
detector  D  when  taking  account  for  the  relative  time  delay,  2  h/c,  where  h  is  the 
pathlength  and  c  is  the  speed  of  light,  can  be  written  as 


IdQ1)  —  ( 


(  2  h\ 

K±u(t)  +  K2  u  +  —  1 


>. 


(ID 


where  K\  and  AT  are  real  numbers  detennined  by  the  losses  in  the  two  paths  and  u(/)  is 
the  analytic  signal  representation  of  the  light  emitted  by  the  source  [7]. 


2.5  Chapter  Summary 

In  summary,  the  ISAL  system  may  be  simplified  into  a  one-directional,  GEO- 
object-to-observatory-point-source  propagation.  Atmospheric  turbulent  layers  will  be 
represented  by  the  FFT-  and  Zernike-based  phase  screens  previously  described,  and  will 
be  temporally  evolved  by  their  individual  means  as  detailed  in  Sec.  5.1.2  and  Ch.  3 
respectively.  The  pulse  lengths  dictate  the  duration  of  light  traveling  through  the 
atmosphere,  which  in  turn  directly  effects  the  time  in  which  the  atmosphere  may  evolve. 
It  is  during  these  brief  pulses  that  attempting  to  combine  the  propagated  beam  with  a 
control  beam  will  effectively  measure  the  theoretical  mixing  efficiency. 
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III.  Temporal  Correlation 


This  chapter  builds  upon  the  understanding  of  Zernike-based  phase  screen 
generation  discussed  in  Sec  2.3.2,  providing  the  ability  for  statistically  accurate  temporal 
evolution  of  a  phase  screen  not  limited  to  the  frozen  flow  approximation.  Atmospheric 
layers  are  modeled  into  the  simulation  by  two  dimensional  phase  screens  representing 
respective  layers  of  turbulence.  This  alternating  free-space  diffraction  and  phase  screen 
accumulation  represents  the  effect  of  propagating  through  an  extended  volume  of 
turbulence.  The  Zemike-based  approach  is  utilized  because  traditional,  FFT-based  phase 
screens  are  limited  to  frozen  flow  evolution.  Without  determining  the  factor  of 
correlation  between  individual  pixels,  the  FFT  methods  can  only  shift,  not  evolve,  with 
statistical  accuracy.  In  the  application  of  IS AL,  the  atmosphere  cannot  shift  any 
significant  distance.  Yet  atmospheric  effects  are  still  observed,  thus,  the  need  to 
investigate  a  new  model  that  permits  turbulent  flow.  This  chapter  develops  the  means  to 
generate  a  phase  screen  that  evolves  as  it  shifts  to  overcome  the  shortfall  of  frozen  flow 
limitations. 

3.1  Expanding  the  Random  Walk  Algorithm 

The  benefit  of  the  Zernike-based  method  is  that  once  a  single  random  screen  is 
generated,  the  correlation  between  each  pixel  is  inherent  within  the  Zemike  polynomials. 
Cain's  [3]  work  with  atmospheric  tilt  has  already  generated  a  statistically  accurate 
algorithm  to  temporally  evolve  tilt,  the  second  and  third  Zemike  polynomials.  By 
expanding  this  algorithm  from  these  two  polynomials  to  any  N  maximum  number  of 
Zernike  polynomials  and  utilizing  the  methodology  of  the  expanded  algorithm,  the 
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generation  of  a  temporal  correlation  of  the  n  vector  is  computed.  Updating  the  n  vector 
is  much  more  computationally  efficient  than  it  is  to  update  the  aj  vector,  as  the  n  vector  is 
independent  and  identically  distributed  (IID).  With  a  temporally  correlated  n  vector,  the 
aj  vector  can  evolve  accurately  and  in  turn,  contiguous  phase  screens  may  evolve.  The 
expansion  of  the  random  walk  algorithm  begins  by  utilizing  Bayes'  theorem  of  the 
conditional  density  function  Eq.  (12), 

f  fvU ,n  _  /nfa+l),n(tfe)(x>y) 

/Wi)i„(«)(*ly)  /n(ti)(y)  <12) 

in  which  the  joint  density  function  f  (n(tk+\),n(tk))(x’  V)  's  zero  mean  Gaussian,  tk  and  tk+1 
are  two  different  observation  times,  and  n(tk )  is  a  single  element  of  the  n  vector  at  time 
tk.  By  denoting  E[n(tk )2]  =  E[n(tk  +  l)2]  as  Rn( 0)  and  E[n(tk)n(tk  +  1)]  as 
Rn(5t),  where  8t  =  tk  +  1  -  tk,  the  joint  density  assumes  the  standard  bivariate 
Gaussian  fonn  [2]: 


fn(tk+l),n(tk)  (x>  y) 


■  x  e 


[  ~Rn(  0) 

y2  2  Rn(St)xy  x2 

^2(R2(0)-R2(<5t)) 

Rn(  0)  S2(0)  Rn(  0) 

(13) 


2tt[1?2(0)  -  R^St)]2 


The  marginal  density  function  fn(tk )  (y)  has  the  standard  fonn  of  a  Gaussian  random 
variable  as  provided  in  Eq.  (14). 


fn(tk)(y ) 


e-y2/2fi„(0) 

■\J  2TcRn  (0) 


(14) 


By  substituting  Eq.  (13)  and  Eq.  (14)  into  Eq.  (12)  the  expanded  conditional  density 
function  is  redefined  as 
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Rn (0) 


fn(tk+l)\n(tk)  (Ay) 


2n(R2(0)  -  R2(8t))/Rn(0) 


x  e' 


>(/?£(0)-/?£(St)) 


(y_vMM\2 

\  Rn( 0)  ) 


(15) 


From  this  expanded  conditional  density,  the  conditional  mean,  Eq.  (16),  and  the 
conditional  variance,  Eq.  (17),  of  n(tk  +  1)  given  n(t/c)  necessary  to  update  the  n  vector 


are  extracted: 


_n(tk)Rn(8t ) 

Mn(tfc+1)|n(tfc)  —  R  (0)  ^ 

2  _(/tf(0) -/?*($,))  (17) 

an(tn+1)\n(tk )  Rn(0) 

Of  significance  at  this  point  is  the  mathematical  representation  of  the  correlation  of  the  n 
vector  in  time  which  once  computed  provides  the  ability  to  update  the  phase 

screens  while  maintaining  both  correlation  and  randomness.  The  conditional  mean  in  Eq. 
(16)  can  be  rewritten  as  an  expectation:  Mn(tfc+1)|n(tfc)  —  £’[n(t2)|n(ti)].  By  substitution 
of  the  inverse  definition  of  the  n  vector  from  Eq.  (6),  n  =  <E>_1a;-,  the  correlation  between 
elements  of  the  n  vectors  can  therefore  be  defined  as 


RnX5t)  =  E 


N  \  /  N 

^  J  (  ^  gy.f1  cfi(t2) 

k= 1  /  \  1=1 


(18) 


Solving  for  the  accurate  correlation  requires  substitution  of  the  a’ s  defined  by  [3]  which 
expands  Eq.  (18)  by  the  definitions  of  oA/c  and  ahi  in  Eq.  (19)  and  Eq.  (20),  respectively: 

/  /  4(w,s)<p;/c(w,  s)0atm(w,s)dwds 


0/,i(t2)  = 


/ J  A(w,  s)(pj  k2(w,s)dwds 

f  f  A(w,s)(pjil(w,s)0atm{w  +  (vx)St,s  +  (vy)St)dw2ds2 
f  f  A(w,  s)q)j  k2(w,  s)dwds 


(19) 


(20) 
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The  aperture  is  represented  as  A  and  the  random  atmospheric  phase  screen  as  0atm,  each 
with  the  same  coordinate  system  (w,s).  Note  that  a/  fe  is  considered  to  be  at  t  =  0  and  a/i( 
is  at  time  t  =  0  +  8t,  where  8t  is  the  difference  in  time.  Finally,  vx  and  vy  are  wind 
velocities.  With  this  substitution,  observatory  parameters,  and  atmospheric  conditions 
are  included  in  the  definition  of  Rn(St)  provided  in  Eq.  (21): 


RnA8t)  =  E 


Is* 


\fc= 1 


f  f  ACw^sJtp/cCw,  s)0atm(w1,s1)dw1ds1 
f  f  A(w,s)<pk2(w,  s)dwds 


J  f  A(w2,s2)<pi(w2,s2)0atm(w2  +  vxSt,s2  +  vySt)dw2ds2 


(21) 


f  f  A(w,  s')(pl2(w,  s)dwds 

The  previous  notation  of  the  indexed  Zernike  polynomials  Zj  has  been  redefined  as  (pk 

and  cpi  in  Eq.  (21),  pennitting  the  update  operation,  which  is  the  means  for  temporally 

correlated  evolution.  By  rearranging,  recognizing  that  the  expectation  operation  is  linear, 

and  grouping  similar  terms,  Eq.  (21)  is  redefined  as  Eq.  (22): 

Vy  _i  — !  J  J  J  J  A(w1,s1)A(w2,s2')(pj:k(w1,s1)(pj,i(w2,s2') 

Uj  1  —  ,l  j  j  A(w,s)q>j  k2(w,s)dwds  f  f  A(w,s)(pj  2(w,s)dwds  (22) 

x  T[0atm(w1,s1)0atm(w2  +  vxSt,s2  +  vxSt)]dw1ds1dw2ds2 

This  is  possible  due  to  the  first  line  in  Eq.  (21)  being  deterministic.  Recognizing  that  the 

expectation  in  the  second  line  is  a  correlation,  Eq.  (22)  is  redefined  as  Eq.  (23)  by  a 

change  of  variables,  tx  —  ( w2  —  wq^Ty  —  (s2  —  %)  and  dzx  —  dw2,  dry  —  ds2. 

-1  5  f  f  5  A(w1,S1)A(w1  +  T X,S1  +  T y)<P;,|c(w1,S1)<P;j!(w1  +  TX,S1  +  Ty) 


RnJ(5t)  =  YJY,=Jk~1=Jl 


J  J  A(w,  s)q>jk2(w,  s)dwds  f  f  A(w,s)<pj  i2(w,s)dwds 


(23) 


X  Ro  [(TX  +  (vx  +  V„)8t,  Ty  +  (yx  +  v^Sjjdw^s^dTy 


The  first  line  of  the  Eq.  (23)  numerator  is  simplified  as  two  independent  correlations  of 
Zemike  polynomials  with  the  aperture,  Eq.  (24)  and  Eq.  (25): 
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(24) 


G  =  A(w1  +  T  X,S1  +  Ty)(pjil{w1  +  T  X,S±  +  Ty)  ( 

Using  Fourier  transform  correlation,  this  is  subsequently  redefined  as  the  inner  double 
integral  operation  P^f. 


Pk,l{jx'Xy)  ~  J  J  A(w±,  S^cpj  k(w±,  S1)A(w1  +  TX,S1  +  T y)<Pj,l(w1  +  T X,S1  +  Ty)  dw1ds1  (26) 

CO  CO 

Pu,  provided  in  Eq.  (26),  is  deterministic  with  the  previously  mentioned  inputs,  and  is 
therefore  only  calculated  once  while  computing  a  By  defining  Pkj  in  this  manner, 

the  quadruple  integral  in  Eq.  (23)  is  simplified  in  Eq.  (27)  as  the  correlation  operation  is 
independent  of  the  outer  double  integral: 


Pnj(St)  = 


JV 

II 

k= 1  1=1 


0.  .  10.  , 


J  J fif.;(r*,Ty)/?g[( Tx  +  vxSt>xy  +  vxSt)]d Txdty 
J  f  A(w,s)(pk2(w,s)dwds  j  J  A(w,s)(pl2(w,s)dwds 


(27) 


Unfortunately  there  is  no  method  to  compute  the  correlation  between  atmospheric  phase 


screens  Rg  directly  using  Kolmogorov  models  for  the  atmosphere  because  Rg( 0)  ->  oo. 
Therefore,  the  relationship  of  the  phase  structure  to  phase  correlation  defined  by  Eq.  (28), 
is  substituted  in  [3]: 


£n,(A)  =  2R^  (0)  -  2R^ («5t).  (28) 

Mathematical  manipulation  is  required  to  use  this  relationship;  it  is  necessary  to  rearrange 
Eq.  (27)  as  well  as  add  zero  (+  [0,0]  -  ffe[0,0]).  This  substitution  redefines  Eq.  (27)  as 


Eq.  (29),  now  in  terms  of  D- (5t)  instead  of  (St): 


N  N 

F>ni8t)= yy®  ^ 

"  fr1  -m  -a 


-1®  -1 


ffp  (  s  I  Ro[o.o]-Rs[o,o]  +  Ro[rx.xy\ 

JJUiU,.uj2  D\  (  +  *  ^  ]drJr'’ 


(29) 


J  J  A(w,  s)^>Jt2(w,  s)dwds  J  J  A(w,  s^cp^fw,  s)dwds 
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By  pairing  the  phase  variances  with  the  appropriate  phase  correlations,  the  correlation 
functions  in  Eq.  (29)  can  be  expressed  in  terms  of  a  phase  structure  function  difference: 


iV 


_!  J  J  Pw(T*’Ty)K(T*  +  vxSt’Ty  +  vxSt )  -  De(jx'Ty)\dTxdTy 


k= 1  /=! 


*  =.  W,  s)<Pk2(  w,  s)3was  f  Ja(  W,  s)<P2  {w,  s)dwds 


(30) 


With  Eq.  (30)  in  terms  of  a  difference  in  phase  structure  function,  Eq.  (3 1)  providing  the 


methodology  to  analytically  calculate  De,  Dn-  may  now  be  solved  for,  and  subsequently 
Rfl,  by  rearranging  the  relationship  of  Eq.  (28),  Rf  =  (D^  —  2)/2.  Recall  that  the  n 
vector  is  zero  mean  and  unit  variance  and  as  such,  f?n(0)  =  1.  Note  that 


Do(jx>1y)  =  6'88 


hi  +Ty\516 
V  ro  ) 


(31) 


3.2  Temporal  Evolution  of  Zernike-based  Phase  Screens 

With  Rf[  calculated  for  a  desired  number  of  polynomials  N,  and  inputs  A,  ro,  and 
wind  speeds  vx  &  vy,  the  random  n  vector  can  be  accurately  updated  with  conditional 
mean,  Eq.  (16)  and  conditional  variance,  Eq.  (17).  The  Rf  vector  only  needs  to  be 
generated  once  for  each  set  of  conditions.  Subsequently,  from  the  definition  of  the 
Zernike  polynomial  coefficient  amplitudes,  Eq.  (6),  the  phase  screen  <p(Sp,  6)  defined  in 
Sec  2.3.2  can  be  temporally  evolved  very  quickly.  Note  that  an  independent  /?„■  is 
required  for  each  individual  phase  screen  given  various  wind  velocities.  The 
computational  requirements  to  generate  a  Rf  vector  is  a  function  of  the  number  of 
Zernike  polynomials  used  to  simulate  a  phase  screen.  For  instance,  it  requires  less  than  a 
minute  to  generate  Rf  with  128  polynomials  on  the  computer  system  used  to  perfonn  the 
analysis  in  this  thesis.  The  system  requires  roughly  12  minutes  to  generate  a  R a  with 
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1,024  polynomials,  the  amount  used  in  all  subsequent  analysis.  As  alluded  to  previously, 
the  tradeoff  between  this  Zernike-based  approach  and  the  FFT-based  approach  is  time. 

To  generate  the  first  8  screens  takes  just  over  an  hour  due  to  the  eight  computations  of 
Rfl.  However,  once  the  /?n-  vectors  are  generated,  the  Zemike-based  simulations 
described  in  subsequent  chapters  are  completed  in  well  under  an  hour.  Inversely,  the 
FFT-based  simulations  generate  the  initial  8  large  phase  screens  in  less  than  a  minute,  but 
the  temporal  simulations  used  for  comparative  analysis  requires  multiple  hours  to 
complete.  This  is  in  large  part  due  to  the  computational  requirements  necessary  to 
perfonn  the  sub-pixel  shifts  of  the  frozen  flow  model  based  on  the  Fourier  shift  theorem. 
A  detailed  description  of  implementing  this  sub  pixel  shift  is  provided  in  Sec.  5.1.2. 

3.3  Chapter  Summary 

While  computationally  it  is  inefficient  to  regenerate  the  R %  vector  each  time  a 
simulation  is  performed,  generating  and  storing  a  database  of  vectors  for  given  input 
conditions  permits  the  efficient  temporal  evolution  of  phase  screens  with  highly  accurate 
correlation.  In  contrast  to  the  frozen-flow  method  in  which  very  large  phase  arrays  must 
be  generated,  stored,  and  then  shifted,  this  Zernike-based  turbulent  flow  method  requires 
only  the  sparsely  populated  correlation  matrix  C,  the  correlation  vector  R^,  and  the 
Zernike  polynomials  Z,  to  be  stored.  The  update  then  becomes  an  extremely  fast  matrix 
operation  which  does  not  require  the  enormous  amount  of  memory  and  processing  time 
required  to  store  and  shift  giant  phase  arrays. 
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IV.  Model  Evaluation 


This  chapter’s  focus  is  the  assessment  of  the  new  Zernike-based  temporal  model 
developed  in  the  previous  chapter  through  laboratory  experiment.  The  first  section 
describes  the  laboratory  configuration  within  which  the  new  model  is  assessed  as  well  as 
the  approach  of  extracting  the  observed  r0  in  the  laboratory.  Next,  the  comparative 
methods  of  evaluating  recorded,  simulated,  and  theoretical  data  are  evaluated.  Finally, 
the  results  of  this  comparison  are  presented. 

4.1  Laboratory  Experiment  Configuration 

The  evaluation  of  this  new  turbulent  flow  model  is  critical  to  the  analysis 
accuracy  of  the  thesis  objectives.  Evaluation  requires  the  comparison  of  recorded 
experimental  data  to  generated  simulation.  The  experimental  data  was  recorded  by 
capturing  one  hundred  successive  photos,  14  photos  per  second,  of  a  pinhole  point  source 
propagated  though  a  concentrated  turbulent  medium  directly  in  front  of  a  2  mm  limiting 
aperture  and  then  focused  by  a  lens  of  focal  length,  //  =100  mm.  A  large  temperature 
differential  generated  the  turbulent  medium.  Extracting  the  rO  of  the  turbulence  in  this 
experiment  is  the  subject  of  Sec.  4.2. 


turbulent  medium 


Figure  5:  Experimental  configuration. 
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4.1.1  Point  Source 


The  point  source  is  created  by  puncturing  a  piece  of  aluminum  foil  with  a  22 
gauge  wire.  The  point  source  was  illuminated  by  a  LED  with  wavelength  A  =  550  mn. 
The  backlit  point  source  was  given  uniform  lighting  to  prevent  directional  propagation  by 
placing  partially  translucent  material  directly  on  the  aluminum  foil.  The  LED  light 
source  was  then  enclosed  to  prevent  ambient  light  noise. 

4.1.2  Optics 

Nyquist  sampling  requires  a  sampling  frequency  to  be  greater  than  or  equal  to 
twice  the  spatial  frequency  cut  off,/,  >2  fc.  The  spatial  frequency  cut  off  for  circular 
apertures  is  equal  to  the  aperture  diameter  A  divided  by  the  wavelength  times  the  focal 
length  f,  given  as  fc  =  A/Af  [6].  Therefore  Nyquist  sampling  is  fs>  lAD.f.  By 
rearranging  and  accounting  for  the  pixel  size  (A)  as  the  spatial  cut  off,  the  limiting 
aperture  diameter  A  can  be  used  to  force  Nyquist  sampling,  as  provided  in  Eq.  (32): 


A  < 


Kfi ) 

2(d) 


(550nm)(0.1042m) 
2(16  pm) 


<  2.32 mm. 


(32) 


The  focusing  length  was  set  based  on  a  97”  propagation  to  the  lens  which  has  a  focal 
length  of  100  mm,  or  3.937”.  Therefore,  by  the  lens  maker’s  equation, 


1  =(  1  M 

3.937"  V4.103"  97"/ 


(33) 


the  focusing  length^  was  set  to  the  image  distance,  approximately  4.1”,  previously 
provided  as  0.1042m. 


The  camera’s  CCD  converts  photons  to  charge,  which  is  then  read  out  serially 
through  register  operations.  The  camera  attempts  to  remove  the  Poisson  distributed 


21 


incident  bias  from  ambient  noise  by  subtracting  the  average  background  from  the  entire 
signal.  See  Fig.  6  for  a  physical  picture  of  the  Camera,  limiting  aperture,  and  focusing 
lens  configuration. 


Figure  6:  Camera  and  optics. 


4.1.3  Turbulent  Medium 

Had  a  turbulence  wheel  been  utilized,  this  analysis  would  be  incorrectly 
evaluating  a  frozen  flow  approximation.  Therefore  a  realistic  concentrated  turbulent 
medium  was  produced  by  using  a  soldering  iron  placed  immediately  in  front  of  and 
below  the  limiting  aperture.  The  iron  has  a  rating  corresponding  to  a  tip  temperature  of 
250  °C.  This  temperature  differential  creates  a  very  strong  turbulent  medium  with 
convective  air  flow.  The  velocity  in  the  center  of  the  convection  vc  of  the  flow  is  given  in 
Eq.  (34)  and  is  a  function  of  the  acceleration  of  gravity  g,  vertical  distance  from  the 
surface  /,  ambient  room  temperature  te,  surface  temperature  ts,  and  the  temperature 
differential  between  the  two  dt=  ts  -  te  [4] : 
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(34) 


vc 


gxlxte  ' 
(273  +dt). 


1/2 


Based  on  a  230  °C  temperature  difference  over  a  5  mm  vertical  distance,  the  vc  is 
calculated  to  be  2.68  cm/sec,  which  translates  to  just  under  2  mm  per  captured  image. 
This  is  reflected  in  the  subsequent  results  section  of  this  chapter  with  small  yet  tangible 
correlations  in  time.  Based  upon  the  long  exposure  OTF  analysis  provided  in  Sec.  4.3  it 
appears  that  the  induced  turbulent  medium  is  Kolmogorov. 


4.2  Extracting  ro 

To  accurately  simulate  the  new  model,  the  first  concern  is  to  match  the  laboratory 
configuration  as  closely  as  possible.  Simulating  this  experiment  is  relatively  easy  given 
the  aperture  size,  vertical  air  velocity,  propagation  distances,  etc.,  except  that  r0  is 
initially  unknown.  Knowing  ro  is  crucial  for  accurately  comparing  simulation  to 
measured  data,  therefore  two  independent  approaches  are  utilized  in  its  extraction. 

4.2.1  Maximum  Correlation 

The  first  method  measures  the  point  spread  function  (PSF)  of  the  raw  data  which 
is  compared  to  theoretical  PSFs  as  a  function  of  r0  through  a  maximum  correlation 
estimator.  The  mathematical  representation  for  cross-correlation  of  real-valued  signals 
across  N  discrete  shifts  is: 


Pit) 


1  y1  (d(0  ~  d)(j(t  +  t)  -  y) 


±y 

N 


r= 0 


(35) 


Here,  d(t )  is  the  measured  intensity  value  and  y(t  +  r)  is  the  value  at  time  t  of  the 
reference  waveform  shifted  by  time  r,  with  their  standard  deviations,  a ^  and  Gy,  to 
nonnalize  the  value  [5].  The  cross  correlation  value  p  reaches  an  absolute  maximum 
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when  the  two  waveforms  are  at  an  optimized  alignment.  Fig.  9  in  Sec.  4.2.3  shows  the 
resulting  plot. 

4.2.2  Gerchberg-Saxton  Phase  Retrieval 

Additionally,  as  a  verification  of  the  maximum  correlation  estimator,  the  recorded 
long  exposure  optical  transfer  function  (OTF)  is  extracted  using  a  Gerchberg-Saxton 
phase  retrieval  algorithm,  and  compared  to  theoretical  long  exposure  OTFs  as  a  function 
of  rn.  The  Gerchberg-Saxton  algorithm  iteratively  uses  the  forward  and  inverse  Fourier 
transforms  to  approximate  the  amplitudes  at  the  source  and  target  plane,  and  the  intensity 
of  the  target  to  the  source  intensity.  This  iterative  process  progressively  estimates,  or 
extracts,  the  phase  as  it  cycles.  Figure  7  depicts  the  entire  cycle. 


PSF  Estimate  <  Sqrt(PSF) 

I  H(w,s) 


Phase  Estimate 

;<p(w\s) 


sj  H(w,s) 


V//new(vv,s)e^(w's) 


Field  at  the  Pupil 

Inverse 
Fourier 
transform 


A' (w,s)elvt'w-s) 


Replace  H  with  new 
PSF  Estimate 


Replace  amplitude.  A, 
with  imaging  system  pupil 


A(w,s)elv("'s'> 

Fourier 

transform 

Figure  7:  Gerchberg-Saxton  algorithm. 

The  square-root  of  the  PSF  estimate  ^Jh(w,  s )  is  combined  with  the  phase 
estimate  e-^O.s)  l0  jn jtiatc  the  algorithm.  This  iterative  process  updates  the  phase 
estimate  as  it  cycles.  Figure  8  shows  the  extracted  phase  from  this  process,  both  wrapped 
and,  by  a  least  squares  error  algorithm,  unwrapped.  The  resulting  mean  squared  error 
(MSE)  plot  between  the  measured  and  theoretical  long  exposure  OTFs  is  presented  in 
Sec.  4.2.3  Fig.  10. 
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Extracted  phase 


LS  Unwrapped  Phase 


Figure  8:  Phase  extracted  from  raw  data. 


4.2.3  Results  of  Extracting  r0 

Figures  9  and  10  have  maximum  and  minimums  at  3.4  mm  and  3.3  mm, 
respectively.  With  these  two  independent  approaches  of  measuring  ro,  it  is  with 
confidence  to  infer  that  the  laboratory  seeing  parameter  is  accurately  estimated  as  3.4 
mm.  This  extracted  ro  allows  for  simulation  of  the  frozen-  and  turbulent-flow  models  for 
comparative  analysis.  Laboratory  comparison  next  compares  the  measured  results  with 
generated  and  theoretical  results. 


Figure  9:  Maximum  correlation  r0  fitting. 
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Figure  10:  Mean  squared  error  r0  fitting. 


Note  that  after  reaching  the  maximum  correlation  and  minimum  MSE  values,  these 
graphs  become  asymptotic  because  as  the  respective  ro  increases,  the  two  comparison's 
approaches  become  diffraction  limited. 

4.3  Comparing  Recorded,  Simulated,  and  Theoretical  Tilt  Correlation  in  Time 

With  a  measured  ro ,  [3]  provides  the  methodology  to  analytically  calculate  the 
theoretical  tilt  correlation  given  input  parameters  A=2  mm,  ro=3.4  mm,  and  wind 
velocities  vx  =  0  &  vy  =  0.00192  m/s.  Based  on  the  concentration  of  the  turbulence 
analysis,  the  atmospheric  phase  is  accepted  as  a  thin  lens  approximation.  The  theoretical 
basis  for  calculating  x  and y  tilt  in  time  is  provided  in  Eq.  (36)  [3]: 


(36) 


Tilt  ax  y  is  calculated  from  the  second  and  third  Zemike  polynomials  cpX:y  which 
correspond  to  x  and  y  tilt  respectively,  the  aperture  A,  and  instant  draws  of  successive 
Zemike-based  turbulent  flow  and  FFT-based  frozen  flow  phase  screens  6atm.  Since  tilt  is 
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zero  mean  [3],  to  calculate  the  temporal  correlation  of  these  tilts,  it  is  important  to 
recognize  that  the  correlation  of  tilt  in  time  is 

i<  (-2)  E [^x,y (Ti)>  @x,y C^-2 )] ■  (37) 

By  Monte  Carlo,  an  average  of  100  simulations  produce  a  Zemike-  and  Fourier-based 

generated  mean  tilt  correlation  in  time.  These  simulations  draw  on  the  same  low  wind- 
velocity  conditions  calculated  in  Sec.  4.1.3  and  the  laboratory  configuration  of  Sec.  4.1. 
The  simulated  measurements  are  compared  to  both  recorded  and  theoretical  values,  all  of 
which  are  presented  in  Table  1.  The  results  indicate  that  the  turbulent  flow  Zernike- 
based  approach  closely  imitates  both  the  recorded,  and  theoretical  tilt  correlation  in  time. 
Inversely,  the  frozen  flow  results  poorly  reflect  both  theory  and  recorded  measurements. 


X  Tilt  Correlation  in  Time 

Y  Tilt  Correlation  in  Time 

Zernike 

Recorded 

Theory 

Fourier 

Zemike 

Recorded 

Theory 

Fourier 

0.1546 

0.1596 

0.1489 

0.8194 

0.1089 

0.1022 

0.1080 

0.0512 

Table  1 :  Tilt  correlation  in  time 


This  table  is  divided  left  to  right  by  the  observed  x  and  y  tilt  correlations  in  time.  The  x 
direction  is  horizontal,  and  the  y  direction  is  vertical  while  viewing  from  the  direction  of 
propagation.  Each  section  is  further  divided  by  the:  Zemike- based  turbulent  flow 
simulation,  laboratory  measured  and  Recorded,  analytically  calculated  Theory  values,  and 
Fourier- based  frozen  flow  modeled  tilt  correlations  in  time. 

In  addition  to  these  findings  on  temporal  tilt  correlation  in  time,  including  1 ,024 
Zernike  polynomials  provides  comparatively  equivalent  high  special  frequency  content  as 
that  of  Fourier-based  approaches.  This  is  evidenced  by  the  simulated  long  exposure  OTF 
of  this  new  turbulent  flow  model  being  equivalent  to  the  OTF  of  standard  frozen  flow 
models  [11].  This  related  technical  paper  shows  that  with  using  a  mere  400  Zernike 
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polynomials  acceptably  mimics  the  frequency  content  of  the  long  exposure  OTF  as 
evidenced  by  Fig.  1 1  which  provides  MSE  values  of  1.14e-4  and  9.08e-5  between  the 
theoretical  and  the  Fourier  and  Zernike  OTFs,  respectively. 


Measured  &  Theoretical  OTFs  Measured  &  Theoretical  OTFs 


Figure  11:  Theoretical  OTF  compared  to  Fourier  (left)  &  Zernike  (right). 


4.4  Chapter  Summary 

In  summary,  comparison  of  the  Zernike-based  temporally  evolving  screen  is 
accomplished  by  pairing  a  simulated  thin  lens  approximation  with  the  parameters  as 
described  in  Sec.  4.1.  The  laboratory  configuration  circumvents  the  need  to  account  for 
multiple  layers  of  atmospheric  turbulence  and  high  wind  speeds,  which  frozen  flow 
approximations  already  accurately  model.  Through  this  laboratory  experiment  it  is 
shown  that  the  new  turbulent  flow  Zemike-based  model  provides  better  tilt  correlation  in 
time  than  the  frozen  flow  approximation  does.  This  finding  is  the  result  of  Monte  Carlo 
analysis  which  compares  theoretical,  recorded,  and  the  two  approaches  to  simulate  values 
of  tilt  correlation.  Further  testing  is  required  to  verify  higher  order  aberrations  are 
correctly  simulated  since  the  laboratory  environment  was  limited  to  an  A  /  ro  <  1 . 
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V.  Simulation  Analysis  and  Results 


This  chapter  discusses  the  methodology  for  accomplishing  the  objectives  outlined 
in  Sec.  1.2.  Based  upon  simulated  phase  screen  analysis,  this  chapter  discusses 
atmospheric  attributes  and  related  approximations,  as  well  as  how  both  the  Zemike-  and 
FFT-based  phase  screens  are  appropriately  simulated.  Additionally,  the  approach  for 
calculating  modulation  efficiency  is  included.  Finally,  this  chapter  provides  the 
analytical  results  for  both  approaches  of  atmospheric  impacts  on  long  pulse  LADAR 
beams  in  the  form  of  mixing  efficiencies. 

5.1  Simulation  Environment 

In  order  to  analyze  what  effects  the  atmosphere  has  upon  mixing  efficiency  as  the 
LADAR  pulse  length  changes,  propagation  through  atmospheric  turbulence  must  be 
accurately  simulated.  The  3.67  m  class  AEOS  telescope  at  the  AMOS  on  Mt.  Haleakala 
is  the  intended  observatory  to  utilize  the  ISAL  capability.  Therefore  utilizing  the  basic 
parameters  of  the  AEOS  will  best  address  the  objectives  of  this  thesis.  Note  that  the  AO 
system  is  not  directly  modeled  into  the  simulation  because  the  AO  cannot  account  for  the 
temporal  effects  as  previously  discussed.  However,  it  is  assumed  that  the  AO  system 
effectively  corrects  wavefront  error  at  time  t  =  0,  or  at  the  initiation  of  each  illumination 
pulse. 

Two  additional  and  forefront  simplifications  of  this  simulation  are  the  absence  of 
anisoplanatic  effects  and  scintillation  noise.  Anisoplanatic  effects  are  of  primary  concern 
while  viewing  objects  large  enough  for  light  from  the  target  to  propagate  through 
different  parts  of  the  atmosphere  to  the  observatory.  For  example,  looking  at  the  moon 
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versus  looking  at  a  distant  star.  Due  to  the  significant  distance  and  a  relative  point  source 
target  anisoplanatic  effects  may  be  ignored  for  this  analysis,  which  greatly  simplifies  the 
simulation.  Secondly,  ignoring  scintillation  noise  is  a  conscious  calculation.  The 
concentration  of  this  analysis  is  to  simulate  and  analyze  idealistic  conditions  so  that 
research  efforts  will  better  understand  mixing  efficiencies.  By  ignoring  scintillation 
effects,  this  analysis  provides  a  baseline  for  such  understanding. 

5.1.1  Atmosphere 

Atmospheric  turbulence  is  represented  by  phase  screens  generated  using  two 
methods:  the  turbulent  flow  model  described  in  Ch.  Ill  and  then  compared  in  Ch.  IV,  and 
the  frozen  flow  model  generated  by  Sec.  2.3.1  and  evolved  by  Sec.  5.1.2.  Eight 
independent  screens  are  generated  representing  thin  lens  approximations  of  the  eight 
respective  atmospheric  regions  shown  in  Fig.  12. 


Figure  12:  Atmosphere’s  structure. 
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These  phase  screen  representations  of  atmospheric  regions  are  the:  Thennosphere, 
Mesopause,  Mesosphere,  Stratopause,  Stratosphere,  Tropopause,  and  Troposphere,  and 
an  eighth  screen  is  dedicated  to  the  fast  flowing,  narrow  current  of  the  Jet  Stream;  which 
is  located  near  the  transition  between  the  Troposphere  and  Tropopause  [1].  The  FFT- 
based  frozen  flow  approach  is  simulated  for  comparative  analysis  to  the  Zemike -based 
turbulent  flow  approach,  and  the  results  of  this  experiment  are  presented  in  Sec.  5.2.2. 

As  a  consequence  of  ignoring  anisoplanatism,  the  atmospheric  phase  screens  do 
not  require  screen  to  screen  propagation.  As  such,  the  eight  independently  evolving 
Zernike-based  screens  depicted  in  Fig.  13  and  FFT-based  screens  shown  in  Fig.  14  are 
summed  together,  the  results  of  which  are  presented  in  Fig.  15  and  Fig.  16,  respectively. 
By  this  summation  a  single  far  field  propagation  is  easily  perfonned  for  each  time  sample 
effectively  modeling  turbulence  as  being  directly  in  front  of  the  telescope.  Due  to  the  40 
Mm  propagation  the  Fraunhofer  diffraction  condition  in  Eq.  (38)  is  easily  satisfied,  i.e., 

z  »  k —  +V  ^max  =  4  x  107m  »  5.712  x  106  m  (38) 

2 

where  z  is  the  40Mm  distance  propagated,  k  is  the  wave  number,  lull  where  X  is  550  nm, 
and  i^J))  are  the  Cartesian  coordinates  in  the  parallel  plane  in  the  x  and  y  directions.  The 
variables  <f  and  rf  denote  the  distance  in  the  x  and  y  directions.  With  this  condition  met, 
a  simple  two-dimensional  FFT  of  the  GEO  target  may  be  perfonned,  accurately  modeling 
a  Fraunhofer  propagation  [7]. 
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Thermosphere 


Mesopause 


Figure  13:  Zernike  phase  screens  representing  atmospheric  layers. 
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Figure  14:  FFT  phase  screens  representing  atmospheric  layers. 
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Figure  15:  Total  Zernike-based  phase,  representing  8  atmospheric  layers. 


Figure  16:  Total  FFT-based  phase,  representing  8  atmospheric  layers. 
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5.1.2  Sub  Pixel  Shifting 

As  previously  stated,  FFT-based  screens  are  temporally  evolved  by  the  frozen 
flow  approximation.  However  due  to  the  time  scales  of  ISAL,  winds  in  excess  of  30 
m/sec  translate  to  simulated  screen  shifts  significantly  less  than  a  single  pixel.  If  this 
analysis  were  simplified,  and  these  fractional  shifts  simply  rounded,  the  screens  would 
not  move.  Therefore,  to  accomplish  a  non-discrete  shift  in  a  discrete  environment  the 
very  large  FFT  screens  are  shifted  using  the  Fourier  shift  theorem.  This  theorem  states 
that  if  jF{g(x,y)}  =  G(fxfi),  then  with  real-valued  constants  a  and  b 


T{g(x  —  a,y  —  b)}  =  G(fx,  fY)e-^xa^b\ 


(39) 


In  essence,  the  shift  in  the  spatial  domain  is  completed  by  introducing  a  linear  phase  shift 
in  the  frequency  domain.  This  linear  phase  shift  e~j2n^xa,^Yb^  is  calculated  by  Euler’s 
formula.  For  instance,  a  13  m/sec  x-directional  wind  velocity  translates  to  a  sub  pixel 
shift  according  to 


X, 


shift 


(13  m/sec)(At ) 

c W: j 


=  0.0145, 


(40) 


where  At  is  1  //sec  (the  time  between  samples),  A  is  3.67  m  (the  primary  aperture 
diameter),  and  N  is  4,096  (the  number  of  pixels  on  the  phase  screen  in  the  x  direction). 
This  fraction  is  fed  into  Euler’s  fonnula 


eJXshift  =  (cos  xshift  +  j  sin  xshift ),  (41) 

which  is  then  factored  into  th efx  and  fy  of  the  Fourier  shift  theorem  for  x  and  y  directional 
winds  respectively.  This  ensures  that  the  spatial  domain  of  the  original  screen  remains 
unaltered,  but  simply  shifted. 


35 


5.1.3  Coherent  Illumination  Beam 


The  ISAL  beam  will  be  utilizing  a  chirped  waveform;  however  this  is  modeled 
into  the  simulation  environment  with  a  rectangle  function  waveform.  The  GEO  space 
object  is  modeled  as  a  single  pixel  in  the  source  plane,  illuminated  by  this  unit  amplitude 
rectangle  waveform  of  varying  length.  This  analysis  pennits  the  ISAL  beam  length  to 
vary  from  10  //sec  to  100  //sec,  by  1  //sec  intervals,  for  a  total  of  91  different  pulse 
lengths.  The  modulation  efficiency  is  calculated  by  the  total  intensity  received  per  pulse, 
Eq.  (46). 

5.1.4  System  Integration 

The  simulation  is  built  into  Matlab®  in  the  following  manner.  First  the  constant 
variables  such  as  aperture  diameter  Fried's  constant  tq,  sampling  sizes  for  the  aperture 
and  receiver  planes  are  set,  and  the  first  1,024  Zernike  polynomials  are  loaded.  From  the 
Zernike  polynomials  and  the  methodology  detailed  in  Ch.  II,  the  covariance  matrix  C  is 
generated.  Random  wind  velocities  for  each  atmospheric  layer  are  randomly  populated 
from  a  uniform  distribution  of  -35  m/s  to  +35  m/s  wind  speeds  in  the  x  and  y  directions, 
permitting  360°  of  directional  wind  of  randomly  varying  intensity  for  each  independent 
atmospheric  layer. 

These  parameters  are  then  fed  into  the  temporal  correlation  algorithm  developed 
in  Ch.  Ill,  providing  the  correlation  vector  Rfi(.5t )  for  evolving  the  eight  Zemike-based 
atmospheric  phase  screens.  Transversely,  these  random  wind  speeds  feed  into  the  frozen 
flow  approach  which  equate  to  the  sub-pixel  shifts  detailed  in  Sec.  5.1.2.  The  Zernike 
phase  screens  are  then  randomly  generated  by  Eq.  (5)  and  Eq.  (6),  the  FFT  screens  are 
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generated  by  the  approach  described  in  Sec.  2.3.1,  and  then  both  are  evolved  in  their 
respective  manner. 

These  eight  independently  generated,  evolving  and  flowing  screens  are  then 
summed  over  the  third  dimension  for  each  time  interval,  where  8t  —  1  [isec  since 
anisoplanatic  effects  are  ignored.  A  single  FFT  propagation  is  perfonned  to  the  aperture 
through  that  instance  of  the  total  atmosphere.  Then  a  second  FFT  propagation  is 
performed  from  the  aperture  to  the  receiver  plane,  where  the  test  and  reference  beams  are 
combined  as  discussed  in  Sec.  2.4  and  are  shown  in  Fig.  17  and  Fig.  18.  The  raw  images 
on  the  left  of  these  figures  are  the  total  test  beam  image  in  the  receiver  plane;  whereas, 
the  combined  images  on  the  right  are  the  test  and  reference  beams  combined.  These 
beams  are  used  to  measure  the  mixing  efficiencies  as  calculated  in  Sec  5.2.  It  is  visually 
apparent  that  this  method  of  heterodyne  detection  greatly  increases  the  SNR. 


Raw  Image  Combined  Image 


50  100  150  200  250  50  100  150  200  250 

Figure  17:  Atmospheric  effects  over  10  juste  pulse. 
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Raw  Image  Combined  Image 


50  100  150  200  250  50  100  150  200  250 


Figure  18:  Atmospheric  effects  over  100  jusec  pulse. 


5.2  Modulation  Efficiency  as  the  Pulse  Length  Varies 

Building  on  the  introduction  of  modulation  efficiency  provided  in  Sec.  2.4,  the 
total  intensity  /  in  both  beams  is  provided  in  Eq.  (42)  and  Eq.  (43),  of  which  gcontroi  and 
gtest  are  the  resulting  fields  in  the  receiver  plane  with  coordinates  (x,v): 


N  N 


l 


control 


-II  1 9  control^  i  9  control^  >y^)\ 


x=l y= 1 


Itest  Ix=lIy=l\9test(.X>  4”  9  control^  >y^\  ■ 


(42) 


(43) 


The  gbeam  s  are  calculated  by  Eq.  (44)  and  Eq.  (45),  the  FFT  of  the  image  in  the  aperture 
plane  Pa,  multiplied  by  the  aperture  A,  and  the  complex  phase  at  the  aperture  plane: 

drest  =  ?{Pa  x^x  ejGat^}  (44) 

(45) 


9 Control  ~P{Pa  x  A  x  j 


Note  that  for  the  control  beam,  atmospheric  perturbations  0atmA,  is  a  zero  matrix  because 
it  is  assumed  that  the  control  beam  used  to  boost  the  SNR  ratio  does  not  incur  any  phase 
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delays  or  aberrations  within  its  delay  loop.  The  test  beam  datmz i  is  the  total  atmospheric 
change  for  a  given  pulse  length  as  previously  described.  The  ratio  of  these  total  energies, 
provided  in  Eq.  (46),  provides  the  modulation  efficiency  M : 


M  = 


*(test) 


(j ’(control)/ '2) 


0(control)/2) 


(46) 


The  approach  to  this  ratio  isolates  the  modulation  efficiency  from  the  total  intensity. 

5.2.1  Total  Atmospheric  Change 

The  total  atmospheric  change  datmA  varies  as  the  pulse  length  varies.  For  shorter 
pulse  lengths,  the  atmosphere  has  less  time  to  change,  thus  the  impacts  the  atmosphere 
has  on  any  light  propagating  through  it  will  understandably  be  less.  Of  particular  interest 
is  the  nature  in  which  the  two  approaches  vary  as  they  evolve  during  the  pulse  length  as 
shown  in  Fig.  19.  The  image  on  the  top  is  the  Zernike-based  approach  whereas  the 
bottom  is  the  FFT-based  approach.  Both  images  are  the  simulated  total  atmospheric 
change  during  a  100  //sec  pulse.  Recall  that  it  is  assumed  the  AEOS  AO  system  corrects 
for  the  atmospheric  wavefront  at  time  t  =  0. 
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Total  Atmospheric  A 


Total  Atmospheric  A 


Figure  19:  Total  atmospheric  change  in  Zernike-based  (top)  and  FFT-based  (bottom). 

Of  additional  concern  are  the  significant  differences  in  amplitude  between  the  two 
approaches.  Note  the  z-axis  scale  in  units  of  radians,  see  Figures  20  and  2 1 . 
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Total  Atmopheric  A  during  pulse 


Total  Atmopheric  A  during  pulse 


Figure  20:  Zernike-based  0atmA  over  10  /rsec  pulse  (top)  and  100  /usee  pulse  (bottom). 
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Total  Atmospheric  A  during  pulse 


Total  Atmospheric  A  during  pulse 


Pixels  u  u  Pixels 


Figure  21:  FFT-  based  0atmA  over  10  jusec  pulse  (top)  and  100  /tsec  pulse  (bottom). 
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The  z  axis  scales  from  the  mesh  plots  of  these  two  figures  indicate  that  between  the  two 
approaches,  the  total  atmospheric  change  is  considerably  less  for  the  frozen  flow  method 
than  that  of  the  turbulent  flow.  It  is  expected  that  large  amount  of  high  frequency  content 
in  Fig.  19  and  Fig.  21  is  largely  due  to  the  sub-pixel  shifting  method  utilized. 

5.2.2  Mixing  Efficiency  Results 

The  differences  in  amplitude  of  atmospheric  change  are  confirmed  and  validated 
by  Monte  Carlo  experiment  of  the  simulation  process  of  Sec.  5.1.  Appendix  A  provides 
the  full  set  of  captured  mixing  efficiencies  for  the  Zernike  results,  and  Appendix  B 
provides  the  corresponding  FFT  results.  The  mean  values  of  the  captured  data  are  plotted 
in  Fig.  22  and  represent  the  average  mixing  efficiency  values. 


0  20  40  60  80  100  120 

Pulse  length  (//sec) 

Figure  22:  Average  mixing  efficiency  with  logarithmic  trend  lines. 
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Two  items  are  readily  apparent  from  this  scatter  plot:  the  FFT  based  approach  has  a  much 
higher  standard  deviation  from  its  trend  and  is  on  average  much  higher  than  the  Zemike 
based  results.  The  a2  for  the  FFT  method  is  2.74%  and  0.35%  for  the  Zernike  based 
method.  Additionally,  the  frozen  flow  approximation  mean  is  on  average  47%  higher 
than  what  the  turbulent  flow  model  predictions.  This  indicates  that  if  the  ISAL  program 
limits  its  investigation  of  mixing  efficiency  to  a  frozen  flow  approximation,  they  will 
potentially  be  over  estimating  mixing  efficiencies  by  roughly  50  percent. 

5.3  Chapter  Summary 

This  chapter  addressed  the  approach  to  the  simulation  used  to  meet  the  objectives 
set  in  Sec.  1.2.  It  built  upon  the  theory  and  validation  of  Ch.’s  III  and  IV,  respectively, 
and  showed  that  as  the  pulse  length  extends,  the  modulation  efficiency  is  degraded  for 
both  the  FFT  and  the  Zernike  approaches.  Of  particular  significance,  the  predicted 
average  values  of  the  frozen  flow  approximation  results  are  greater  than  expected.  This 
chapter  showed  that  for  small  time  scales  the  atmosphere  can  in  fact  significantly  evolve, 
reaffirming  the  understanding  that  the  atmosphere  does  not  in  fact  adhere  to  frozen  flow. 
Additionally,  it  was  shown  that  the  new  Zemike  model  provides  more  tightly  fitting 
predictions  given  the  measured  standard  deviation  of  the  mean  mixing  efficiencies 
between  the  two  approaches. 
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VI.  Conclusions  and  Recommendations 


This  chapter  summarizes  the  efforts  of  this  research  and  the  resulting  findings. 
Particular  contributions  are  highlighted  and  the  challenges  are  detailed.  In  addition, 
recommendations  for  future  work  are  included  to  help  further  the  research  of  turbulent 
flow,  temporally  evolving  phase-screen  generation. 

6.1  Summary 

As  discussed  in  Ch.  1,  the  objective  of  this  research  has  been  to  evaluate  the 
atmospheric  effects  on  long  pulse  LADAR  beams.  To  meet  this  objective,  this  thesis 
predominately  focused  on  the  development,  comparison,  and  utilization  of  an  entirely 
new  temporally  evolving  Zemike-based  turbulent  flow  phase  screen  generation  model. 
Building  on  the  understanding  of  generating  independent  Zernike -based  phase  screens 
established  in  Ch.  II,  the  temporal  model  was  developed  by  expanding  a  random  walk 
algorithm.  This  random  walk  algorithm  was  previously  applied  to  atmospheric  tilt  [3], 
and  Ch.  Ill  expanded  the  algorithm  to  include  any  N  number  of  Zemike  polynomials. 

This  model  was  then  validated  through  laboratory  experiment  in  Ch.  IV  and  used  to 
provide  the  analysis  of  the  objectives  listed  in  Sec.  1.2.  The  simulation  application  of  the 
mathematics  behind  the  temporal  aspect  to  Zemike-based  phase  screens  was  by  far  the 
most  difficult  aspect  of  this  research  effort.  The  analysis  of  the  objectives  was 
completed  in  Ch.  V  and  compared  to  FFT-based  frozen  flow  analysis.  It  was  determined 
that  as  the  pulse  length  elongates,  the  mixing  efficiency  declines.  Furthermore,  it  was 
shown  that  the  frozen  flow  approach  has  a  higher  standard  deviation  from  its  trend  line, 
while  the  turbulent  flow  approach  provided  much  tighter  fitting  data  points.  The 
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assumptions  made  at  the  initiation  of  this  research  effort,  such  as  the  limitation  of  the 
frozen  flow  approximation  to  relatively  small  atmospheric  effects,  were  continued. 

6.2  Conclusions  of  Research 

It  was  detennined  that  as  the  pulse  length  extends,  mixing  efficiencies  decline. 
Based  upon  the  Zemike  results,  a  roughly  8%  efficiency  at  10  jj.  sec  and  2%  mixing 
efficiency  at  100  //sec  closely  follows  a  logarithmic  decline.  Similar  in  general  decline, 
but  strikingly  different  in  noise  levels,  the  frozen  flow  approach  potentially  overestimated 
the  mixing  efficiencies  of  the  new  turbulent  flow  model.  Additionally  the  turbulent 
flow  model  was  validated  through  laboratory  experiment.  The  experiment  captured 
images  of  a  point  source  with  a  concentrated  turbulent  medium  directly  in  front  of  the 
limiting  aperture,  and  through  a  PSF  likelihood  maximization  algorithm  and  Gerchberg- 
Saxton  phase  retrieval  algorithm,  extracted  realistic,  turbulent,  temporally  evolving 
phase.  Modeling  the  laboratory  conditions  into  simulation  and  then  applying  both  the 
new  Zemike-based  turbulent  flow  model  and  the  FFT-based  frozen  flow  model,  it  was 
detennined  that  the  turbulent  flow  outperformed  the  frozen  flow  in  accurately  predicting 
tilt  in  time.  Through  previous  and  directly  related  analysis  [1 1],  it  was  shown  that  this 
new  Zemike  approach  outperforms  the  frozen  flow  approximation  in  predicting 
atmospheric  effects  in  temporal  analysis. 

6.3  Significance  of  Research 

As  stated  previously,  the  true  significance  of  this  research  was  the  development  of 
the  new  turbulent  flow  model  of  atmospheric  phase  screen  generation  and  its  comparison. 
This  new  model  outperformed  the  frozen  flow  approximation  in  comparing  both 
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laboratory  results  and  theoretical  tilt  correlations  in  time.  With  these  results,  a  baseline 
may  be  struck  for  future  efforts  within  the  ISAL  program.  Furthermore,  the  new  model 
may  be  applied  to  an  assortment  of  applications  in  Space  Situational  Awareness  (SSA), 
Defense  Advanced  Research  Programs  Agency  (DARPA),  and  deep  space  surveillance, 
to  name  a  few.  The  predominate  application  to  the  relevant  programs  directly  involves 
more  accurately  estimating  temporal  effects  when  imaging  through  the  atmosphere.  The 
novelty  of  this  new  model’s  approach  is  of  considerable  significance,  as  currently  there  is 
no  other  temporal  evolution  phase  screen  generation  model  that  is  not  limited  to  frozen 
flow  approximations,  much  less  any  such  model  that  has  been  validated  as  this  one  has. 

6.4  Recommendations  for  Future  Research 

The  overarching  question  in  response  to  proposing  this  new  approach  to  phase 
screen  generation  has  been,  “Can  this  methodology  be  applied  to  Fourier-based  methods 
of  phase  screen  generation?”  The  answer  is  yes.  However,  the  number  of  calculations 
for  correct  correlation  will  balloon  to  an  N2  number  of  calculations  for  individual  pixel 
updates.  Instead  of  building  on  the  inherent  correlation  between  Zemike’s,  it  is  expected 
that  the  pixel  to  pixel  correlation  update  be  rooted  within  the  power  spectrum  used  for 
generating  the  individual  FFT-based  phase  screens.  Instead  of  generating  an  R„(St) 
vector,  an  R„(S,)  matrix  must  be  generated  for  each  set  of  given  atmospheric  conditions, 
observatory  parameters,  etc.  This  correlation  matrix  would  then  feed  into  the  conditional 
mean  and  conditional  variance  update  matrices.  Finally,  in  similar  fashion  to  the 
considerably  easier  simulation  described  in  the  previous  chapters,  the  correlation  matrix 
can  be  paired  with  FFT  based  phase  screen  generation,  and  each  pixel  could  then  be 
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updated.  The  theoretical  mathematics  to  accomplish  this  is  not  overly  difficult  and  will 
likely  follow  a  very  similar  approach  to  that  of  Ch.  III.  However,  the  difficulty  in 
adapting  such  math  into  simulation  space  is  known  to  be  considerable. 

6.5  Chapter  Summary 

This  chapter  summarized  the  activities  of  this  research  and  the  results  obtained. 
Accomplishing  the  objectives  presented  in  Ch.  I  supported  the  hypothesis  that  as  the 
LADAR  pulse  elongates,  the  atmospheric  effects  increase.  This  degradation  in  mixing 
efficiency  appeared  to  follow  a  logarithmic  decline.  The  challenges  overcome  and 
significant  contributions  were  also  highlighted.  Finally,  this  chapter  outlined  the 
expected  approach  to  developing  a  Fourier-based  temporal  model  not  limited  to  frozen 
flow  approximations. 
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Appendix  A  -  Zernike-based  Mixing  efficiency 
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0.0378  0.0467  0.0291  0.0409  0.0345  0.0355  0.0495  0.0431  0.0328  0.0289 
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0.0363  0.0287  0.0329  0.0151  0.0152  0.0303  0.0145  0.0253  0.023  0.0253 

0.0339  0.0289  0.0267  0.0341  0.0484  0.0124  0.0528  0.0295  0.0305  0.0287 

0.0256  0.0355  0.0285  0.0375  0.0344  0.0304  0.0302  0.0309  0.0274  0.0271 

0.0484  0.0358  0.0457  0.0335  0.0444  0.0405  0.024  0.032  0.0442  0.0372 

0.03702  0.03294  0.03432  0.03098  0.033  0.029  0.02984  0.02856  0.03058  0.02614 

80  nsec  81  nsec  82  nsec  83  nsec  84  nsec  85  nsec  86  nsec  87  nsec  88  nsec  89  nsec 

0.022  0.0295  0.0301  0.0278  0.0216  0.0222  0.0266  0.0337  0.0342  0.0184 

0.0117  0.0323  0.0265  0.0318  0.0229  0.0217  0.0263  0.0277  0.0303  0.0185 

0.039  0.022  0.0227  0.034  0.0412  0.0283  0.0174  0.0209  0.0262  0.0226 

0.0325  0.0148  0.0249  0.0283  0.0264  0.0408  0.0269  0.03  0.0197  0.0208 

0.0357  0.0386  0.0288  0.0254  0.0365  0.0264  0.0363  0.0312  0.0355  0.037 

0.02818  0.02744  0.0266  0.02946  0.02972  0.02788  0.0267  0.0287  0.02918  0.02346 

90  nsec  91  nsec  92  nsec  93  nsec  94  nsec  95  nsec  96  nsec  97  nsec  98  nsec  99  nsec 

0.0118  0.0055  0.0453  0.0197  0.0196  0.0122  0.021  0.0234  0.0173  0.0305 

0.0169  0.0251  0.0294  0.0303  0.0224  0.0207  0.0212  0.0229  0.0134  0.0212 

0.0172  0.0061  0.0168  0.0187  0.0184  0.0337  0.0182  0.0195  0.0181  0.0223 

0.0499  0.0254  0.0265  0.0374  0.0254  0.0258  0.028  0.0285  0.0425  0.0108 

0.0298  0.0354  0.0393  0.0327  0.0368  0.0311  0.0422  0.0232  0.0323  0.0201 

0.02845  0.023  0.028  0.029775  0.02575  0.027825  0.0274  0.023525  0.026575  0.0186 

100  nsec 

0.0069  0.0206  0.0206  0.0226  0.0308  0.02365 
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Appendix  B  -  FFT-based  Mixing  efficiency 


10  nsec 

11  nsec 

12nsec 

13  nsec 

14  nsec 

15  nsec 

16  nsec 

17  nsec 

18  nsec 

19  nsec 

0.1383 

0.106 

0.1068 

0.104 

0.0932 

0.0988 

0.0936 

0.1158 

0.0837 

0.1035 

0.1062 

0.1179 

0.136 

0.093 

0.045 

0.0402 

0.044 

0.0716 

0.09 

0.072 

0.1447 

0.124 

0.078 

0.071 

0.0424 

0.042 

0.0396 

0.0244 

0.016 

0.026 

0.1756 

0.1658 

0.1149 

0.0876 

0.087 

0.082 

0.077 

0.0664 

0.0618 

0.0776 

0.1062 

0.1395 

0.1064 

0.0924 

0.0908 

0.0772 

0.0556 

0.0596 

0.0888 

0.0808 

0.1342 

0.13064 

0.10842 

0.0896 

0.07168 

0.06804 

0.06196 

0.06756 

0.06806 

0.07198 

20  nsec 

21  nsec 

22  nsec 

23  nsec 

24  nsec 

25  nsec 

26  nsec 

27  nsec 

28  nsec 

29  nsec 

0.0972 

0.1242 

0.0922 

0.0535 

0.0621 

0.0702 

0.0605 

0.05265 

0.0580 

0.0999 

0.0904 

0.0963 

0.0796 

0.0787 

0.0666 

0.0661 

0.0670 

0.0621 

0.0715 

0.0706 

0.054 

0.0693 

0.0882 

0.0562 

0.054 

0.0846 

0.1025 

0.0625 

0.063 

0.1030 

0.0711 

0.0891 

0.0895 

0.0661 

0.0607 

0.0711 

0.0715 

0.0576 

0.0666 

0.0841 

0.0324 

0.0345 

0.0455 

0.0310 

0.0895 

0.0918 

0.0846 

0.0877 

0.063 

0.027 

0.06903 

0.08271 

0.07902 

0.05715 

0.0666 

0.07677 

0.07722 

0.06453 

0.06444 

0.07695 

30  nsec 

31  nsec 

32  nsec 

33  nsec 

34  nsec 

35  nsec 

36  nsec 

37  nsec 

38  nsec 

39  nsec 

0.0885 

0.1251 

0.0835 

0.0279 

0.0445 

0.1044 

0.0855 

0.0927 

0.0747 

0.0725 

0.0918 

0.0495 

0.0522 

0.0755 

0.1053 

0.0522 

0.0639 

0.0517 

0.0832 

0.0735 

0.0715 

0.0243 

0.0072 

0.0189 

0.0589 

0.0925 

0.0936 

0.0855 

0.0679 

0.0670 

0.0958 

0.0382 

0.0373 

0.0198 

0.0234 

0.0508 

0.0198 

0.0555 

0.0355 

0.0193 

0.0576 

0.094 

0.0935 

0.0994 

0.0364 

0.0441 

0.0576 

0.0963 

0.1015 

0.0486 

0.08025 

0.06633 

0.05463 

0.04824 

0.05373 

0.06876 

0.06408 

0.07632 

0.07254 

0.05616 

40  nsec 

41  nsec 

42  nsec 

43  nsec 

44  nsec 

45  nsec 

46  nsec 

47  nsec 

48  nsec 

49  nsec 

0.0187 

0.0665 

0.0515 

0.1107 

0.0702 

0.1134 

0.0675 

0.0333 

0.1075 

0.0705 

0.0202 

0.0486 

0.0855 

0.1314 

0.0345 

0.0423 

0.0657 

0.0595 

0.0375 

0.0441 

0.0395 

0.1105 

0.126 

0.1089 

0.1377 

0.0351 

0.0105 

0.0765 

0.09 

0.0666 

0.1005 

0.1062 

0.0702 

0.0335 

0.0485 

0.0765 

0.0955 

0.0531 

0.0585 

0.0345 

0.0513 

0.0475 

0.0655 

0.0285 

0.054 

0.0477 

0.0385 

0.0725 

0.0625 

0.0819 

0.04595 

0.07605 

0.07965 

0.08262 

0.06768 

0.063 

0.05607 

0.0594 

0.07056 

0.06066 

50  nsec 

51  nsec 

52  nsec 

53  nsec 

54  nsec 

55  nsec 

56  nsec 

57  nsec 

58  nsec 

59  nsec 

0.074 

0.08 

0.0196 

0.0196 

0.0364 

0.0464 

0.0744 

0.02 

0.046 

0.0844 

0.0776 

0.0928 

0.0276 

0.0276 

0.0148 

0.0716 

0.0384 

0.0136 

0.0684 

0.0908 

0.0532 

0.0492 

0.0504 

0.0504 

0.034 

0.0792 

0.0196 

0.044 

0.0708 

0.0624 
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0.0684 

0.0544 

0.0348 

0.0348 

0.0492 

0.0088 

0.0428 

0.0428 

0.0408 

0.0464 

0.0504 

0.0644 

0.0972 

0.0972 

0.0384 

0.0848 

0.0696 

0.0516 

0.0412 

0.0728 

0.06472 

0.06816 

0.04592 

0.04592 

0.03456 

0.05816 

0.04896 

0.0344 

0.05344 

0.07136 

60  nsec 

61  nsec 

62  nsec 

63  nsec 

64  nsec 

65  nsec 

66  nsec 

67  nsec 

68  nsec 

69  nsec 

0.0785 

0.0525 

0.098 

0.0365 

0.0441 

0.0476 

0.0567 

0.0155 

0.0399 

0.0385 

0.0615 

0.0651 

0.0595 

0.0715 

0.0505 

0.0505 

0.0555 

0.035 

0.0511 

0.0189 

0.0185 

0.0056 

0.0735 

0.0294 

0.0644 

0.0553 

0.0658 

0.0505 

0.028 

0.0665 

0.0437 

0.0444 

0.0574 

0.0505 

0.0658 

0.0322 

0.0469 

0.0259 

0.0845 

0.0065 

0.0147 

0.0445 

0.0238 

0.0511 

0.0371 

0.0525 

0.0441 

0.0749 

0.0651 

0.0505 

0.04354 

0.04242 

0.06244 

0.04907 

0.05243 

0.049 

0.05383 

0.04018 

0.05369 

0.03619 

70  nsec 

71  nsec 

72  nsec 

73  nsec 

74  nsec 

75  nsec 

76  nsec 

77  nsec 

78  nsec 

79  nsec 

0.0567 

0.0387 

0.0096 

0.0474 

0.0525 

0.0663 

0.0363 

0.0801 

0.0183 

0.021 

0.0531 

0.0726 

0.0663 

0.0582 

0.0336 

0.0393 

0.0261 

0.0684 

0.0471 

0.0171 

0.0048 

0.0096 

0.0159 

0.0585 

0.0231 

0.0351 

0.0555 

0.0303 

0.0477 

0.0525 

0.0183 

0.0591 

0.0618 

0.0168 

0.0618 

0.0375 

0.0726 

0.0651 

0.048 

0.0216 

0.0606 

0.054 

0.0417 

0.0714 

0.0504 

0.0528 

0.0552 

0.0522 

0.0312 

0.0897 

0.0387 

0.0468 

0.03906 

0.05046 

0.04428 

0.0462 

0.04914 

0.05922 

0.03846 

0.04038 

80  nsec 

81  nsec 

82  nsec 

83  nsec 

84  nsec 

85  nsec 

86  nsec 

87  nsec 

88  nsec 

89  nsec 

0.0723 

0.0396 

0.0282 

0.0354 

0.0447 

0.0513 

0.0612 

0.0519 

0.0231 

0.0555 

0.0849 

0.0525 

0.0165 

0.0426 

0.0285 

0.063 

0.0264 

0.0186 

0.0366 

0.0201 

0.0225 

0.0168 

0.0477 

0.057 

0.0426 

0.0615 

0.0468 

0.0429 

0.0144 

0.0396 

0.0318 

0.0741 

0.0537 

0.0567 

0.0414 

0.0564 

0.0459 

0.0522 

0.0564 

0.0345 

0.0342 

0.0258 

0.0537 

0.0339 

0.0483 

0.072 

0.0519 

0.021 

0.0612 

0.0219 

0.04914 

0.04176 

0.03996 

0.04512 

0.0411 

0.06084 

0.04644 

0.03732 

0.03834 

0.03432 

90  nsec 

91  nsec 

92  nsec 

93  nsec 

94  nsec 

95  nsec 

96  nsec 

97  nsec 

98  nsec 

99  nsec 

0.0506 

0.0206 

0.032 

0.0332 

0.0304 

0.052 

0.0436 

0.0258 

0.0226 

0.0242 

0.035 

0.0342 

0.0262 

0.0592 

0.0394 

0.0338 

0.0362 

0.0266 

0.0414 

0.0388 

0.0372 

0.0423 

0.0213 

0.0738 

0.0399 

0.0567 

0.0591 

0.0483 

0.0438 

0.0129 

0.0396 

0.0309 

0.0207 

0.0411 

0.0264 

0.0138 

0.0375 

0.0411 

0.0402 

0.0297 

0.0555 

0.051 

0.0363 

0.0693 

0.0393 

0.0423 

0.0177 

0.0318 

0.0366 

0.0225 

0.04358 

0.0358 

0.0273 

0.05532 

0.03508 

0.03972 

0.03882 

0.03472 

0.03692 

0.02562 

100  nsec 

0.0213 

0.0216 

0.0336 

0.045 

0.06 

0.0363 
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